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Abstract 

We present next-to-next-to-leading order corrections to the inclusive production 
of the Higgs bosons at the CERN Large Hadron Collider (LHC) including finite top 
quark mass effects. Expanding our analytic results for the partonic cross section 
around the soft limit we find agreement with a very recent publication by Harlander 
and Ozeren [1]. 

PACS numbers: 12.38.Bx 14.80.Bn 



1 Introduction 



With the launch of the Large Hadron Collider (LHC) at CERN the particle physics ex- 
periments enter a new energy domain with the hope to discover new phenomena. Exper- 
imental observations are expected to provide hints for the still open questions within the 
Standard Model. Among these questions is the mechanism of the electroweak symmetry 
breaking which in most practical theories provides masses to the particles. The tradi- 
tional implementation of the symmetry breaking implies the existence of a new particle, 
the Higgs boson, which so far has not been detected in high-energy collisions. 

For the intermediate Higgs boson mass, favoured by the results of indirect searches, the 
most important production channel at a hadron collider is gluon fusion, gg — > H, mediated 
by a top quark loop. During the last 20 years enormous efforts have been made to evaluate 
higher order corrections to this process. 

The leading order (LO) result has been presented in Refs. [2-5] and already almost 15 
years ago also the next-to-leading order (NLO) QCD corrections became available [6, 
7]. More recently also the next-to-next-to-leading order (NNLO) corrections have been 
evaluated [8-11]. While the NLO results are exact in the top quark and Higgs boson 
masses (in this context, see also Ref. [12]), the NNLO results rely on the effective theory 
built in the limit of the large top quark mass (see, e.g., Refs. [13, 14] for the three- loop 
corrections to the effective ggH coupling). It is well known that this approximation works 
surprisingly well at NLO, leading to deviations from the exact result that are less than 2% 
for Mh < 2M t [15]. It is one of the aims of the present paper to investigate the validity 
of the large top quark mass approximation at NNLO. 

During the last few years there appeared several improvements over the fixed-order cal- 
culation. Among them is the soft-gluon resummation to next-to-next-to-leading [16] and 
next-to-next-to-next-to-leading [17] logarithmic orders and the identification (and resum- 
mation) of certain n 2 terms [18] which significantly improves the perturbative series. We 
furthermore want to mention Ref. [19] where the gluon-gluon channel has been consid- 
ered in the limit of large center-of-mass energy a/§. Recent numerical predictions of Higgs 
boson production in gluon fusion both at the Tevatron and the LHC are summarized in 
Ref. [20]. 

Very recently the effects of the finite top quark mass on the Higgs boson production in 
hadron colliders have been reported in Ref. [1]. The authors of that reference carried out 
an asymptotic expansion of the corresponding production diagrams and obtained results 
for the first four terms in the p = Mfj/Mf expansion, where the secondary expansion 
around the soft limit (i.e. for x = Mjj/s — > 1) has been performed to a sufficiently high 
order. 

In this paper we present the results of an independent calculation of these finite top mass 
effects, confirming the conclusions of Ref. [1] . Similarly to Ref. [1] , we apply asymptotic 
expansion in 1/M t to the full QCD diagrams and evaluate a few first terms of the series. 
However, our results are not expanded near the soft limit and the x-dependence of the 
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cross section (which is valid below the top pair threshold, as discussed further) is retained. 

The remainder of the paper is organized as follows. In the next Section we discuss the 
cross sections in the individual channels of partonic reactions at NLO and NNLO. We 
consider in particular the behaviour near x — > and describe a method to extrapolate our 
results to the limit of large s. In Section [3] we use the results of Section [2] and numerically 
evaluate the hadronic cross section for the LHC and the Tevatron. Our conclusions are 
presented in Section HI 



2 Partonic cross section 



2.1 Notations and calculation details 

We introduce the following notation for the partonic cross section: 
fr a (a {0) 4 ° s A (1) 4- f as Y A (2) 4- 

with 



and ij denoting one of the possible initial states: gg, qg, qq, qq, or qq', where q and q' 
stand for (different) massless quark flavours!]] At NNLO the Higgs boson in the final state 
may be accompanied by zero, one or two gluons, or by a light quark pair. In general, the 
quantities A^ depend on x and p. Leading order mass dependence is then described by 
the function fo(p, 0) given in Eq. (4) of Ref. [21]. 

Factoring out the exact LO top quark mass dependence as in Eq. (CQ) is a common practice. 
In what follows, by the "infinite top quark mass approximation" we mean that only the 
quantities A)- are evaluated for M t — > oo, but A^o remains exact in M t . 

At the LO only AfJ is different from zero and given by 

Ag = 5(1 -x). (3) 

At the NLO the functions A^, Aqg and A^ are not zero, and at the NNLO one has to 
consider all five contributions: A^, A^, A^| , Agq , and A^,. 

The first results on the p dependence of NNLO cross sections appeared in Refs. [21] 
and [22], where the virtual part of A^ was evaluated to 0(p 4 ) and 0(p 2 ), respectively. 
In this paper we adhere to the notations of Ref. [21]; in particular, we use ai 5 ^ and the 
on-shell top quark mass. 



1 It is understood that ghosts are always considered together with gluons. 
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Figure 1: Sample forward scattering diagrams whose cuts correspond to the LO, NLO 
and NNLO corrections to gg — > H. Dashed, curly and solid lines represent Higgs bosons, 
gluons and top quarks, respectively. 

Since at the LO only the virtual contribution is present, we do not discuss it further. In 
order to evaluate the real corrections to A ^ , we exploit the optical theorem and compute 
the imaginary part of the four-point amplitudes ij — > ij. At NLO and NNLO this requires 
the evaluation of three- and four-loop diagrams, respectively. Some sample diagrams are 
shown in Fig. [TJ Note that only the cuts dissecting the Higgs boson propagator and one 
or two massless lines need to be included. 

We generate the diagrams with the help of QGRAF [23] supplemented by additional scripts 
that eliminate the vanishing graphs. At the next step we apply the asymptotic expansion 
(see, e.g., Ref. [24]) in the limit M f 2 3> s, Mfj, implemented in two independent programs: 
q2e/exp [25,26], and an in-house Perl implementation. 

This procedure factorizes the original triple-scale forward scattering functions into mas- 
sive vacuum integrals (with a single scale M t ) up to three loops, that can be evaluated 
with MATAD [27], and four-point one- and two-loop integrals dependent on s and M#. 
After the reduction performed with the Laporta algorithm [28, 29] we arrive to a set of 
master integrals that have been studied in Ref. [10] (Appendix B). Unfortunately, that 
reference contains a number of misprints. We independently evaluated these integrals by a 
combination of soft expansion^] and differential equation methods (a detailed comparison 
to be provided elsewhere). 

Using these techniques we have been able to compute NNLO contributions through 0(p 2 ) 
for the channels gg and qg, and through 0(p 3 ) for qq, qq and qq'. 

2.2 NLO results 

The integral representation of the NLO corrections to the partonic cross section can be 
found in Refs. [6,7]. In analytic form the p° and p 1 terms (gg channel) were presented 
in Refs. [7] and [30], respectively. In Ref. [1] terms up to 0(p 3 ) have been demonstrated. 

2 We acknowledge help with cross checks of the soft expansion by Robert Harlandcr and Kemal Ozeren. 
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Here it is convenient for us to include terms through 0(p A ) given by 
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Figure 2: Partonic NLO cross sections for the (a) gg, (c) qg and (d) qq channel as functions 
of x for Mn = 130 GeV. The expansion in p — > (dashed lines) is compared with the 
exact result (solid lines). Lines with longer dashes include higher order terms in p. The 
interpolation (see text) is shown as a dotted line. In order to demonstrate the smallness 
of the region below threshold we show in (b) the result for the gg channel using the linear 
x-scale. The vertical line at x = M]j/(AMf) ~ 0.14 indicates the threshold for top quark 
pair production. 



In Fig. [2] functions in the three channels are shown depending on x, evaluated for 
Mh = 130 GeV and M t = 173.1 GeV [31], including successively higher orders in p 
(dashed lines). The exact result is plotted as a solid line. One sees that the leading term 
in p is smooth and demonstrates a reasonably good agreement with the exact curveJl 
However, the higher order terms in p introduce divergences at x — > which are the most 
severe for the qq channel. This signifies the breakdown of the assumption that Mf ^> s 
for large s. Note, however, the decent convergence above the threshold for the top quark 

3 For qq the agreement is not that obvious, however, the hadronic results have the proper order of 
magnitude. 



6 



pair production, i.e., for x > x t h = M^/i^Mf). 

In order to reduce the dependence on unphysical divergences near x = at the NNLO, 
where exact cross sections are not known, we here devise two practical recipes and test 
them against the NLO results. 

For the quark channels no information beyond p expansion is available. Thus, considering 
that these contributions are numerically suppressed and that x t h limits the applicability 
of our asymptotic expansion, we use the following 

Option 1: 

For x > x t h, we use the complete result including all known 0(p n ) corrections, and 
for x < x t h the infinite top mass approximation. 

As will be demonstrated in the following section, this introduces an error in the hadronic 
results that is < 50%, which, if also true at the NNLO, has a sufficiently small effect on 
the hadronic cross section compared to the overall scale uncertainty. 

For the primary production channel, gg —> H, NLO and NNLO asymptotics near x — > 
have been found in Ref. [19]. Also, NLO plots suggest that in this channel there are no 
pronounced threshold effects at x = x t h- Thus, here we can use a more educated 

Option 2: 

Complete 0(p n ) result is matched at some point x m < x t h to a function 3Ci + 
ax (NLO) or — 9C2lnx + b (NNLO), where coefficients C\ and C2 are tabulated 
in Ref. [19] and a, b and the matching point chosen to provide the most 

"natural" smooth behaviour of the function. 

Note that this recipe is different from the procedures suggested in Ref. [19] and Ref. [1]. 
The reason is that our results (as in Eq. (J3])) contain genuine l/x n poles in 0(p n ) con- 
tributions that are not present in the infinite top mass result used in Ref. [19] and are 
masked by the soft expansion of Ref. [1] . 

We found that an x m such that the function and its first derivative match smoothly is 
a good choice at the NLO; at the NNLO, matching at x m = Xth/^ produces reasonable 
results for 110 GeV < Mh < 300 GeV and is consistent with the region of x where 
higher 0(p n ) corrections demonstrate good convergence. By varying the constants and 
interpolating function shapes we have checked that the dependence of the hadronic cross 
section on the exact details of the matching procedure is quite small and that only the 
asymptotics near x — > are important. 

In Figs. [2] (a) and (b) the thus obtained NLO extrapolation in the gg channel is shown as 
a dotted line. One observes very good agreement with the exact curve, and the difference 
in the hadronic cross sections is negligible. 
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2.3 NNLO corrections 



Combining the virtual part of Agjj calculated in Refs. [21,22] with the real contributions 

(2) 

we arrive at the 0(p n ) corrections to the quantities A^- for n = 0,1,2 for gg and qg 
reactions, and n = 0,1,2,3 for the remaining channels. Our results expressed in terms 
of harmonic polylogarithms are quite lengthy and can be obtained from the authors on 
request. The O(p ) terms exactly reproduce the expressions found in Ref. [10]. Expanding 
the higher 0(p n ) corrections in (1 — x) <C 1 we find complete agreement with Ref. [1]. 

In Fig. [3] we present non-singular parts of the functions A^ for ij = gg, qg, qq, qq, qq' 
as functions of x. Here one can observe a behaviour similar to that at the NLO: the 
higher order terms in p develop more severe singularities near x — > 0, however, below the 
threshold the results converge. The dotted curves in Figs. [3]^a) and (b) demonstrate the 
"Option 2" extrapolations described above. The further numerical analysis is based on 
these extrapolations. 



3 Hadronic cross section 

The hadronic cross sections are given by the convolution of the partonic cross section 
dij-^H+x with the corresponding parton distribution functions (PDFs), which is conven- 
tionally written as follows!^ 



k,l(H{g,u,...,b,u,...,b} 





~d£ k i 


[ dr 


hills 


. dr . 



(r, (jl f ) & k i^ H+x (s = rs, p F ) , (7) 



where d£fc//dr is the so-called luminosity function given by 

[T,p F ) = dXi dx 2 fk/p{xi,PF)fl/p{x2,PF)S(T-XiX 2 ) 

Jo Jo 



d£ki 
dr 



For the further discussion we adopt a slightly different parametrization in terms of the 
"natural" parameter x = M% / s and the distinct production channels: 

(a;, p F ) frij^H+x (x, p F ) , (9) 

with the straightforward modifications to the corresponding weights. For example, the 

4 In this paper we concentrate onpp collisions at the LHC peak energy y/s = 14 TeV. The modifications 
for pp collisions at the Tevatron are obvious. 
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Figure 3: Partonic NNLO cross sections for the (a) gg, (c) qg, (d) qq, (e) qq, (f) qq' 
channels functions of x for Mh = 130 GeV. Lines with longer dashes include higher order 
terms in p. In (b) we also show the gg channel in the linear scale. The dotted line in (a) 
and (b) corresponds to the matched result. 



quark-gluon luminosity is defined as 
dC 



'<19 



dx 



x,/i F ) = 2 V / dxi / dx 2 fg/p(x 1 ,fj, F )f q /p(x2,n F ) 

, - t, JO JO 



(10) 



q£{u,...,b,u 

sx 



X\X2 



Ml 



sx 



2 ' 



1000 




Figure 4: Luminosity functions for ij = gg (solid), qg (short dashed), qq (long dashed), 
qq (dash-dotted), and qq' (dotted line). 



For illustration we show in Fig. H]the quantities cLCy/cLr for ij = gg, qg, qq, qq, and qq'. 
One notices that for x — > there is a rapid decay of all luminosity functions which is one 
of the main reasons that at the NLO the heavy top approximation works extremely well. 

We use the parton distribution function (PDF) set MSTW2008 [32] and the a s evolution at 
LO, NLO and NNLO when computing predictions to the cross section at the corresponding 
order. To discuss the numerical effect of our calculation we decompose the prediction of 
the total cross section into its LO, NLO and NNLO contributions: 
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and denote the heavy top quark approximation with an additional subscript oo. 

Let us in a first step discuss the channels involving quarks which are treated using "Op- 
tion 1" as described in Section [2j In Fig. [5] we show the M#-dependence of the NLO 
contribution to the hadronic cross section originating from the quantity A^j (cf. Eq. (TO) 
normalized to the exact result as coded in HIGLU [33] for Mh between 110 GeV and 
300 GeV. For the qg channel one observes that the infinite top quark mass approximation 
provides between 40 and 50% of the exact result. After including the p and p 2 term this 



10 




0.2 




0.4 



0.2 







150 



200 
M H (GeV) 

(a) 



250 



300 







150 



200 
M H (GeV) 

(b) 



250 



300 



Figure 5: Ratio of the NLO hadronic cross section including successive higher orders in 
1/M t (from short to long dashes) normalized to the exact result, (a) qg and (b) qq 

is improved to about 60% for the smaller Higgs boson masses whereas for Mh = 300 GeV 
the heavy top quark mass is practically unchanged. Similarly, for the qq channel we ob- 
serve an improvement by about 10 to 15% for the Higgs boson masses around 140 GeV. 
For Mh = 300 GeV the heavy top expansion is practically equivalent to the exact result. 

The analogous curves at NNLO can be found in Fig. [6] where we normalize the result on 
the infinite top quark mass approximation. In all cases the power-suppressed terms lead 
to an increase of the cross section between 4% and 10% for the quark-gluon and up to 
25% for the quark-anti-quark channel in our range of Higgs boson masses. The very rapid 
convergence is observed for the qq and qq' channels where the contribution beyond the 
1/M t 2 term is practically zero. 

Let us finally turn to the numerically most important contribution, the gg channel treated 
with the matching to s — > oo asymptotics denoted above as the "Option 2" . In Figs. [3(a)- 
(c) we demonstrate the NNLO contribution to the hadronic cross section (cf. Eq. ( fill) ) 
normalized to the infinite top quark mass result. The difference is that in (a), the exact 
LO top quark mass dependence is factored out as in Eq. ([T]), while in (b) the partonic 
cross sections both in numerator and denominator are strictly expanded in p. Finally, 
in (c) we expand Alq in the numerator but keep it exact in the denominator. 

For the fully expanded option (b) one observes for Mh = 300 GeV corrections up to 40% 
originating from the linear p term which further increase to almost 60% after including the 
p 2 term. However, when the exact leading-order top quark mass dependence is factored 
out (case (a)), the corrections amount to at most 8%. Considering the fact that the NNLO 
terms contribute about 10% of the total NNLO cross section we conclude that the top 
quark mass suppressed terms at NNLO alter the prediction by less than 1%. This justifies 
the use of the heavy top mass approximation for the evaluation of the NNLO hadronic 
cross section. The latter conclusion is also obtained from Fig. [3(c). 
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Figure 6: Ratio of the NNLO hadronic cross section including successive higher orders in 
1/M t (from short to long dashes) normalized to the infinite top quark mass result, (a) 
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In Fig. 0(d) we also take into account the exact LO and NLO contribution and again study 
the effect of the different p terms. Similar to the case (c), we leave in the denominator 
the exact LO mass dependence and consider various expansion depths of the NNLO 
contribution in the numerator. This plot can be directly compared to the left panel of 
Fig. 7 in Ref. [1]. Very good agreement is observed; the minor differences can be traced 
back to the different matching procedures. 



4 Conclusions 

In this paper we present the NNLO production cross section of the Standard Model Higgs 
boson including finite top quark mass effects. Our calculation is based on the evaluation 
of the imaginary part of the forward scattering amplitudes which, via the optical theorem, 
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Figure 7: (a), (b) and (c): Ratio of the NNLO hadronic cross section (gg contribution) 
including successive higher orders in 1/M t normalized to the infinite top quark mass 
result. In (a) the exact LO mass dependence is factorized both in the numerator and 
denominator. In (b) numerator and denominator are expanded in p, and in (c) only the 
numerator is expanded, (d) shows the prediction of the gluon-induced inclusive Higgs 
production cross section up to NNLO normalized to the heavy top limit. 



directly leads to the total cross section. We apply the asymptotic expansion in order to 
obtain correction terms suppressed by the heavy top quark mass. 

We observe rapid convergence of the series below the threshold for the production of real 
top quarks, i.e. for s < 4M t 2 . However, the region of small x = Mfj/s demonstrates l/x n 
singularities as a consequence of our expansion procedure. For the numerically dominant 
gluon-gluon channel we match our results to the large s limit, curing thus those artificial 
singularities and obtaining stable predictions for the hadronic cross section. 

The numerical impact of the top quark mass suppressed terms is below approximately 
1% and thus about a factor ten smaller than the uncertainty from scale variation. Let us, 
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however, stress that this result was not obvious a priori. Our calculation justifies the use 
of the heavy top quark mass approximation when evaluating the NNLO cross section. 

In addition, we confirm the results of Ref. [10] for the infinite top quark mass, and the 
Mi-suppressed terms calculated in Ref. [1]. 
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